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Abstract 

We perform an all-orders resummation of the QCD Adler D-function for the 
vector correlator, in which the portion of perturbative coefficients contain- 
ing the leading power of 6, the first beta-function coefficient, is resummed. 
To avoid a renormalization scale dependence when we match the resum- 
mation to the exactly known next-to-leading order (NLO), and next-NLO 
(NNLO) results, we employ the Complete Renormalization Group Improve- 
ment (CORGI) approach in which all RG-predictable ultra-violet logarithms 
are resummed to all-orders, removing all dependence on the renormaliza- 
tion scale. We can also obtain fixed-order CORGI results. Including suit- 
able weight-functions we can numerically integrate these results for the D- 
function in the complex energy plane to obtain so-called "contour-improved" 
results for the ratio R and its tau decay analogue Rr. We use the difference 
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between the all-orders and fixed-order (NNLO) results to estimate the un- 
certainty in as{M'^) extracted from measurements, and find as{M^) = 
0.120±0.002. We also estimate the corresponding uncertainty in a{Mz) aris- 
ing from hadronic corrections by considering the uncertainty in R{s), in the 
low-energy region, and compare with other estimates. Analogous resumma- 
tions are also given for the scalar correlator. As an adjunct to these studies 
we show how fixed-order contour-improved results can be obtained analyti- 
cally in closed form at the two-loop level in terms of the Lambert 1^-function 
and hypergeometric functions. 
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1 Introduction 



The correlator of two vector currents in the Euchdean region is a fundamental 
ingredient in constructing a number of inclusive hadronic QCD observables 
of great importance in testing the theory. By taking a logarithmic energy 
derivative one can define the so-called Adler D-function, D{s). By analyti- 
cal continuation to the Minkowski region this quantity can then be directly 
related to the ratio ,R{s), of the total e"'"e~ hadronic cross section to the 
point leptonic cross- sect ion, and also to the analogous ratio Rr of the total 
hadronic decay width of the r lepton normalized to the leptonic decay width. 
The analytical continuation can be elegantly formulated as a contour integra- 
tion of D[s) together with a weight function around a circle in the complex 
energy s-plane |l|, |^. Performing this integration numerically with D{s) ap- 
proximated at some fixed order of perturbation theory then automatically 
resums to all-orders an infinite subset of potentially large analytical continu- 
ation terms involving powers of vr^ and beta-fuction coefficients, which arise 
in the running of the coupling around the integration contour. These terms 
are usually truncated in the direct fixed-order perturbative expansion of the 
Minkowskian quantity. Such approximations are referred to as "contour- 
improved" . 

The remaining uncertainty in these "contour-improved" predictions for R 
and Rr comes from the uncalculated higher order terms in the perturbation 
series for D{s), which has presently only been computed exactly to ©(a^^) in 
the limit of massless quarks |^, and with some approximations including the 
contribution of top and bottom quarks 0] . We shall assume massless quarks 
in these investigations. There are also effects due to non-perturbative power 
corrections 0, |^, |^. We shall focus in this work on the former perturba- 
tive uncertainties. There are two interrelated aspects to these. Fixed order 
perturbation theory predictions have a dependence on the renormalization 
scheme (RS) chosen to define the coupling. In particular they depend on a 
dimensionful renormalization scale fi. Further, given a choice of RS there is 
an uncertainty due to the unknown 0(af) and higher uncalculated pertur- 
bative terms. To estimate this one needs to perform a necessarily approx- 
imate all-orders resummation of these terms. A well-motivated framework 
to accomplish this is provided by the so-called "leading-6" approximation 

1^ (sometimes also referred to as "naive nonabelianization" [9-11]), which 
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amounts to resumming to all-orders the portion of perturbative coefficients 
containing the highest power ofb= 1(11 - 2Nf) , the ffist beta-function 
coefficient for SU(A^) QCD with Nf active massless quark flavours. This 
can be accomplished since in the large- A^j limit one has an exact all-orders 
result for the Adler D-function [12-14]. The leading-6 resummation is then 
performed by replacing Nf by ~ 36). Whilst the leading-6 approach 

is motivated by the structure of renormalon singularities in the Borel plane 
[0, |10|, [Tl[, and also by a QCD skeleton expansion [jl5[, it is effectively just 



the ffist "one chain" term in the skeleton expansion, and does not include the 
multiple exchanges of renormalon chains needed to build the full asymptotic 
behaviour of the perturbative coefficients, and there are no ffim guarantees 
as to its accuracy. The strongest statement that can be made follows from 
an analysis of the operators which build the leading ultraviolet renormalon 
singularity |TB, |T^. One can prove that in the case of the vector Adler func- 



tion the re-expansion of the leading-6 result in powers of Nf correctly gives 
the asymptotics of the portion of perturbative coefficients proportional to 
Nf'^~^N^' in ?T,**^-order perturbation theory, with accuracy 0{l/n) |T^. In 



practice for the exact NLO and NNLO coefficients of the Adler D-function 
the level of agreement of these individual coefficients is at the ten percent 
level, much better than would be expected from the above weak asymptotic 
result, remarkably in the Nf = or large- A^ limit agreement is at the few 
percent level Q. We may therefore hope that the leading-6 approximation is 
indicative of the size of uncalculated higher-order corrections. A remaining 



difficulty, ffist emphasised in |T8[ , is the scale dependence of the leading-6 re- 
summations, if one tries to match them to the known exact NLO and NNLO 
perturbative coefficients. This matching ambiguity means that any result 



may be obtained by varying the scale . It was pointed out that claims f^, [TT 



that comparison of fixed-order and "leading-6" resummed results indicated 
rather large uncalculated perturbative corrections for were undermined 
by this matching problem. In Ref . [|T^, ^ the difficulty was resolved by per- 



forming a leading-6 resummation for the Effective Charge beta-function pO| 
corresponding to D{s). This scheme- invariant construct then unambiguouly 
determines Rr- This approach revealed a rather small uncertainty due to 
uncalculated higher-order corrections. 

In this paper we wish to formulate the resummations in a very closely re- 
lated, but technically much more straightforward way. The renormalization 
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scale, /i, dependence of fixed-order QCD perturbation theory is an artefact of 
the way renormalization group (RG) improvement is customarily performed. 
The two crucial features are the use of a scale /i proportional to the physical 
energy scale, Q, of the process, and the truncation of the perturbation se- 
ries at fixed-order. As argued recently |2l]] one should rather keep /z strictly 
independent of Q. Fixed-order perturbation theory with fi constant does 
not then satisfy asymptotic freedom, and one is forced to sum to all-orders 
the RG-predictable unphysical logarithms of /i and physical ultraviolet (UV) 
logarithms of Q from which the perturbative coefficients are built. This 
so-called "Complete RG-improvement" (CORGI) |21| serves to cancel all /x- 
dependence between the unphysical renormalised coupling and the 

unphysical logarithms of /i in the coefficients, and one directly trades un- 
physical /i dependence for the physical Q-dependence. The idea can also be 
generalized to processes, such as structure function moments, which involve 



a factorization scale as well as a renormalization scale 12111. The CORGI 



approach as formulated in Ref.[^ is exactly equivalent at NLO to the Effec- 
tive Charge approach of Grunberg in which the UV logarithms are also 
completely resummed in exactly the same way, whilst the remaining RG- 
predictable effects are parametrized in a different, but a priori equally rea- 
sonable way. Our plan is to perform a leading-6 resummation for the Adler-Z? 
function in the CORGI approach. As we shall see this is extremely straight- 
forward to implement and the resulting resummed result can be written as 
a sum over exponential integral functions, representing the contributions of 
the ultraviolet and infra-red renormalons in the Borel plane. In contrast the 



resummation of the Effective Charge beta-function in Refs.0, |T9l involved 
a complicated numerical inversion of a function. We anticipate that the two 
approaches should yield very similar results. We shall obtain leading-6 re- 
summed and contour-improved CORGI results for the quantities R{s), and 



Rt-. As in Refs. ||18| , |19| we shall use the difference in fixed-order and resummed 
results to estimate the uncertainty in q;s(M|) obtained from measurements of 
Rr, using more recent experimental data We shall also fit the leading-6 
resummation to the spectral distribution for hadronic r decay [^, We 



shall attempt to estimate the uncertainty in the hadronic corrections to the 
value of the QED coupling at the Z pole, a(M|), using the resummed and 
fixed-order results for R{s) in the energy ranges 5 < y/s < oo GeV, and 
2.8 < ^/s < 3.74 GeV, and using inclusive data in the remaining ranges, as 



in Ref.|2^. We shall compare our result for q;(M|) with that of Ref.[^ 
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which uses standard fixed-order perturbation theory. Using recent large- A^/ 
results on the scalar correlator we shall also perform a contour-improved 
leading-6 resummation for the Higgs decay width. Finally, we show how the 
analytical continuation to the Minkowski region can be performed in closed 
analytical form at the two-loop level in terms of the Lambert M^-function 
and hypergeometric functions. 



The organization of the paper is as follows. In Section 2 we shall review 
the contour integral representation of R{s) and Rr in terms of D{s), and 
describe a simple numerical algorithm for evaluating it. In Section 3 we 
shall review the fixed-order perturbative, and all-orders large- iV/ results for 
D{s) , and show how the leading-6 all-orders resummation in the CORGI 
approach can be written in closed form as a sum of exponential integral 
functions. In Section 4 we estimate the uncertainty in as(M|) obtained from 
Rt measurements, and also fit the leading-6 resummed results to the spectral 
function. In Section 5 we shall estimate the uncertainty in the hadronic 
corrections to a(M|) as discussed above. In Section 6 we perform a contour- 
improved leading-6 resummation for the Higgs decay width, and in Section 
7 discuss how the analytical continuation to the Minkowski region can be 
performed analytically at the two-loop level in terms of the Lambert W- 
function and hypergeometric functions. Section 8 contains a discussion and 
our conclusions. 



2 Contour integral representation of Minkowski 
observables 

We shall mainly be concerned in this work with two inclusive QCD observ- 
ables. The first is the e~^e~ -R-ratio, defined by 

cr(e+e~ hadrons) , . 

In SU(A^) QCD perturbation theory 

R{s) = (l + IcfR{s)^ + (y^Qi^ Ri') ' (2) 
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with Qf denoting the quark charges, summed over the flavours accessible 
at a given energy. Cp is the SU(iV) Casimir Cp = {N^ — 1)/2N, and s 
denotes the timelike Minkowski squared momentum transfer. R denotes the 
perturbative corections to the parton model result. It has the perturbative 
expansion 

R{s) = a{l + 5]r„a") , (3) 

n>0 

where a=as{^'^)/'Tc denotes the RG-improved coupling. R denotes so-called 
"light-by-light" contributions which enter at O(a^). We shall ignore this term 
in our all-orders resummations. The ratio Rr is defined analogously as a ratio 
of the total r hadronic decay width to its leptonic decay width, 

T(t ^ Ut- + hadrons) , ^ 

Rr = — ^ . 4 



Its perturbative expansion has the form 



Rr = N{\V^,f+\Vusf)SEW 



5 aimr"^) ~ 

1 + T77 ' +Rr + 5pc 



(5) 



with Vud and Vus CKM mixing matrix elements. Since the energy scale 
s = lies below the threshold for charmed hadron production only three 
flavours 'U,(i,s, are active. The a(m^) term denotes the leading QED electro- 



magnetic corrections, and S'£;vf~1-0194 ||26[ represents further electroweak 
corrections. 5pc denotes possible power corrections. Rt has a perturbative 
expansion 

R, = a(l + ^r>") . (6) 

n>0 

In this case there are no "light-by-light" corrections because summing over 
u,d and s quarks (Z]Q/)^=0. Both R and Rr can be directly expressed in 
terms of the transverse part of the correlator of two vector currents in the 
Euclidean region, 

{q,.qu - g^uQ^Ms) = ^ttH J rf^xe*^-" < o|r[j^(x)j,(x)]|o > , (7) 

where s = —q^ > 0. In fact it is convenient to take a logarithmic derivative 
with respect to s and define the Adler Z)-function, 

D{s) = -^^n(s) . (8) 
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This can be represented by an expression analogous to Eq.(2) involving D 
and i), where D[s) has the perturbative expansion 

b{s) = a(l + J^rf^a") . (9) 

n>0 

A generic Minkowskian observable R{sq) can then be related to D{—s) by 
analytical continuation from Euclidean to Minkowski. This can be elegantly 
formulated as an integration around a circular contour in the complex energy 



squared s-plane pJ[] , 

R{so) = Wie)D{soe'')d9 , (10) 

ZnJ-TT 

where W{6) is a weight function which depends on the observable R. For 
W{e) = 1 one has R{so) = R{so), and for W{e) = (1 + 2e*^ - 26^^^ - e^*^) 
one has R{m'^)=Rr. If one expands D{soe'^^) as a perturbation series in 
a=a(soe*^) and numerically performs the 6 integration term-by-term one ob- 
tains "contour-improved" perturbative results in which at each order an in- 
finite subset of analytical continuation terms present in the conventional 
perturbation series of Eqs.(3),(6) are resummed. These terms are potentially 
large and involve powers of vr^ and beta-function coefficients, as is easily seen 
by expanding a in powers of a(so) and integrating. In this paper we shall 



focus on this "contour-improved" version of perturbation theory. In Ref. |19 
detailed comparisons of the performance of the two versions were made, and 
the importance of resumming the analytical continuation terms was empha- 
sised. 



An obvious numerical algorithm for evaluating the integral in Eq.(lO) is 
to split the range from 6=0,71 into K steps of size A6 = ir/K and perform a 
sum over the integrand evaluated at ^„=nA^, n=0,l,...,K. So that 

/\0 ^ K ^ 

R(so)-—[W{0)D{so) + 2ReY,W{er,)D{sn)] (11) 

n=i 

where s„=Soe*"^^. In practice we perform a Simpson's Rule evaluation. Writ- 
ing the perturbation expansion for D{sn) we have 

Disn) = dn + dial + didl + ... , (12) 
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where we have defined an=a(s„). An efficient strategy is to start with 
Go = a(so) and use Taylor's theorem step-by-step to evolve a„ to Un+i, using 

a„+i = an- i—bB{a„) — 6^5(a„)5'(a„) + i—— 6^[5(a„)5'(a„) 



2 ' ' 8 ^ ^ 48 



+5(a„)^5"(aO] + 0{M^) + ... (13) 
where B{x) is the truncated beta-function 

B{x) = x^ + cx^ + C2a;^ + . . . , (14) 

so that a satisfies 

da h h 
^ = --(a^ + ca^ + caa^ + ...) = --B{a) . (15) 

Here b = (33 — 2Nf)/6, and c = (153 — 19Nf)/12b are the first two univer- 
sal beta-function coefficients for SU(3) QCD with Nf active massless quark 
fiavours. The higher coefficients Ci,i > 1 are scheme-dependent. The above 
use of Taylor's theorem is much faster to implement than the standard ap- 
proach of solving the integrated beta-function equation with complex 
renormalization scale Sn to find dn at each step. 



3 Fixed-order and resummed expressions for 
D{s) in the CORGI approach 

In the CORGI approach one avoids renormalization scale /^-dependence by 
performing a complete resummation of the ultraviolet logarithms which build 
the dependence of the observable on the physical energy scale This is 
equivalent to directly relating the observable to the dimensional transmuta- 
tion parameter of the theory ||2^, Ajjg say. In this way one can write the 



CORGI series for D{s), 

D{s) = ao{s) + X2al{s) + Xga^ + . . . + X„a;j+i + . . . . (16) 

Here ao(s) is the CORGI coupling which may be written in terms of the 
Lambert W^-function defined imphcitly by W {z)exp(W (z)) = z [EH] as, 



c\\^W{z{s))\ 



9 



where K£,=e'^/^{2c/h)~'^^^ K-jjg^ with d the NLO perturbative coefficient di for 
D{s) in Eq.(9), in the MS scheme with /i^ = s. ao(s) is the couphng in the 
scheme with /x^ = e~'^'^^^s, and the non-universal beta-function coefficients, 
Cj, (i > 1) all zero. In this scheme rfi = 0, and it is exactly equivalent 
at NLO to the Effective Charge approach of Grunberg Standard RG- 



improvement in this scheme completely resums all ultraviolet logarithms, and 
is equivalent to the CORGI approach which can be formulated in any scheme 
X2 is the NNLO scheme-invariant combination 

X2 = C2 + d2- cdi - d\ , (18) 

built from the perturbative coefficients di and ^2 and beta-function coef- 
ficients. The NLO and NNLO coefficients di and ^2 are known exactly 
and so NNLO contour-improved CORGI predictions can be straightforwardly 
obtained for Minkowski observables R{so), using the numerical integration 
described in Section 2. Since ao(s) is known in closed form in terms of the 
Lambert ly-function, which has a well-defined branch structure in the com- 
plex plane, one can evaluate it directly, avoiding the Taylor's theorem trick 
in Eq.(13). In fact one needs to use the W-i branch of the function (in 
the nomenclature of Ref.|2^) on the range of integration [0,7r], and the Wi 
branch on the range [— 7r,0]. As we shall discuss in Section 7 one can, in 
fact, avoid using the numerical Simpson's Rule integration all together for 
the case of the e~^e~ i?-ratio where W{6) = 1, and perform the integral in 
closed form in terms of logarithms of the l^-function. 

In order to assess the likely accuracy of the fixed-order perturbative ap- 
proximation we can attempt to approximate at the present uncalculated co- 
efficients di, {i > 2) in D{s) using the so-called "leading-6" approximation. 
A given coefficient dn can be written as an expansion in powers of Nf, so 
that we have 

d^ = rfWiv; + S^~'^NJ~' + . . . + . (19) 

The large-Nf coefficient d]^^ can be computed exactly to all-orders since it 
derives from a restricted set of diagrams in which a chain of n fermion bubbles 
is inserted in the initiating quark loop [0, |13|. Motivated by the structure 
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of renormalon singularities in the Borel plane one can convert this expansion 
into the so-called leading-6 expansion ||, by substituting Nf = (33/2 — 36), 
to obtain 

rf„ = rfWfe" + + . . . + . (20) 

The leading-6 term d^^=d^^b'^ is then used to approximate dn- Since d^^ = 
(— 3)"(i|^l6" it is known to all-orders. Using the exact large-A'^/- result one 
finds that the explicit all-orders result for d!f^^ in the so-called V-scheme, i.e. 
MS with scale /i^ = e~^^^s, is given by [|l^ 

dL'\V) = ^n!^^[-2n-^ 

+^ E s(l-2-2^)(l-2^— 2C2.+i)]6V (21) 

^ + ■'-f+l>s>0 

The resulting leading-6 resummation 

00 

£)(^) = a(l + E4'^V), (22) 

fc=0 

may then be defined as a principal value (PV) regulated Borel Sum, 

/•oo _ 

D^^\l/a) = PV dze-'/^BlD^^^z) . (23) 

^ 

Here B[D^^^](z), denotes the Borel transform which contains an infinite set of 
single and double poles at z = zi = ^ corresponding to infra-red renormalons 
, IRi, and an infinite set of ultra-violet renormalons, UVi, aX z = —zi. The 
structure is 

The residues at these poles can be computed from the exact all-orders large- 
Nf result. The UV and IR renormalon contributions can then be easily 
expressed in terms of the exponential integral function. 



00 



Ei{x) = - I dt^ , (25) 
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where for IR renormalons x > and one defines Ei{x) by taking the Cauchy 
principal value of the integral. The arbitrariness in regulating the IR renor- 
malon contributions reflects the fact that the perturbation series needs to 
be combined with the power corrections of the operator product expansion 
(OPE) to obtain a well-defined result The absence of a relevant opera- 
tor of dimension two in the OPE for the vector correlator is in accord with 
the fact that the singularity IRi is not present, and the nearest singularity 
to the origin in the Borel plane is in fact UVi, which generates the leading 
asymptotic behaviour p|, 

4^)(y)^il^!i±M^!('_iyV . (26) 



27 V 2 

One can then write the UV renormalon and IR renormalon contributions as 
infinite sums over the Ei functions, 

oo 

+{AoU) - ^jMM} , (27) 

and 

D^(F)|,^ = e-^''z2Bo{2)Et{Fz2) 

oo 

Y.z,{e~^^^EiiFz,)[Fz,iBoij) + z,B^{j)) - z,B,{j)] 

j=3 

-(BoU) + z,B,{j))} . (28) 

Here we have defined F=l/av, where ay is the coupling in the V-scheme. 
The Ao{j),Ai{j) are related to the residues of the UVj poles, with p 

8(-lOM + 2) 4K-l)'--(2, + 3) 



Because of the conformal symmetry |^ of the vector correlator the UV 
residues are directly related to the IR residues with Bo{j) = —Ao{—j) and 
BiU) = 

-Aii-j) for j>2, and Bo{l) = B^{1) = 5i(2) = 0, and Bo{2) = 1 §. To 
evaluate the contour integral in the complex s-plane using this D^^\F) result 
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one needs to modify the definition of the Ei functions to cope with the fact 
that their argument involves l/ay(soe*^) which is complex for nonzero 6. The 
appropriate generalization uses the function Ei{n, z) defined by 



oo 



Ei{n,z) = dt^. (30) 

This function is analytic everywhere in the cut complex z-plane, but has a 
branch cut along the negative real axis. One needs to replace Ei{—Fzj) in 
the UV contribution by —Ei{l, Fzj), and Ei{Fzj) in the IR contribution by 
—Ei{l, —Fzj) + i'Ksign{Im{Fzj)), where the discontinuity across the branch 
cut is removed by the final iir contribution ||^. The final result for D{F) is 
simply the sum of the UV and IR contributions. The sums in Eqs.(27),(28) 
are rapidly convergent since the A{j) and B{j) coefficients have a de- 
pendence for large j. For the numerical results to be reported in Section 
4 we used Nuv = 15 and Njji = 17 terms respectively in the two sums. It 
is sensible to arrange that Njr = Njjv + 2 , since the symmetry properties 
above mean that Ao(j) = —Bo{j + 2), this ensures that the 0(a) term in the 
perturbation series of Eq.(9) has the correct unit coefficient Bq{2) = 1. 

The final step is to use the above results to perform an all-orders re- 
summation in the CORGI approach. We would like to formally perform the 
resummation 

DcoBGi = ao + X2al + . (31) 

n>2 

SO that the exactly known NNLO X2 coefficient is included , with the re- 
maining unknown coefficients approximated by X^^\ X^^\ . . ., the leading-6 
approximations. Note that qq is the full CORGI coupling defined in Eq.(17), 
so that all the RG-predictable ultraviolet logarithms involving the exact NLO 
coefficient di are completely resummed. This resummation is most easily 
achieved by noting that the combination 



po = bin (^j^ - di{n) , (32) 

is scheme-independent. At the leading-6 level the coupling a^^\s) is defined 
by the one-loop formula 

a(^ns) = • (33) 

^ ^ 61n(v^/A) ^ ^ 
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In the CORGI scheme in leading-6 approximation d!{'' = 0, and so by eval- 
uating the invariant combination po in Eq.(32) in the V scheme and the 
CORGI scheme one can relate the couplings in the two schemes, 

^ ^ +dnV). (34) 



AL) (L) 



It then follows straightforwardly that the formal resummation in Eq.(31) is 
given by 

DcoRGi = D^"-^ {- + df\v)] + {X, - Vo , (35) 



in which the D^^'^ term is the all-orders sum with the exact X2 replaced by 
and the second term corrects for this. One can obtain approximate 
N^LO and higher CORGI results by truncating the sum in Eq.(31). The 
X!^^ can be readily calculated by using the leading-6 relation between the 
IZ-scheme and CORGI couplings in (34). One easily finds 



n+l 



(36) 



where the symbol C„[/(a)] denotes the coefficient of a" in the power series 
expansion of /(a). The SnK^) ^e directly generated using the explicit 
result in Eq.(21). 

Using the above results we can now straightforwardly generate all-orders re- 
summed and fixed-order contour-improved CORGI results for the Minkowski 
observables Rr and R. We shall perform some phenomenological studies in 
the next two sections. 



4 Resummed versus fixed-order predictions for 

The ratio Rr defined by Eq.(4) has been the subject of a wide-ranging exper- 
imental study by the ALEPH collaboration |22[. If events involving strange 
quarks are removed from the data, they find Rr = 3.492±0.016. Setting 
Vus = 0, and Vud = 0.9754±0.0007, and estimating the power correction 
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Order of pert. theory 



Figure 1: Fixed-order CORGI results for i?^ in N"LO perturbation theory 
(starred points), compared to the all-orders resummation (solid line) fitted 
to ALEPH data. 



contribution to be 5pc = — 0.003±0.004 [Q, one finds from Eq.(5) the exper- 
imental value Rr = 0.2032l!]:[!J^^. The QED contribution has been neglected. 
One can then obtain all-orders leading-6 resummed and fixed-order contour- 
improved CORGI results as described in Sections 2 and 3. We use Nf = 3 

(3) 

and fix A^;^ so that the all-orders result reproduces the measured central 
value Rr = 0.203. The results are shown in Figure 1. The solid line is the 
all-orders resummed result fixed to the data, and the starred points show the 
N'^LO fixed-order CORGI results. We see that the NNLO (n=2) fixed-order 
result, which is the highest order exactly known, is in rather good agreement 
with the all-orders resummation. The leading-6 approximated N^LO results 
show an oscillatory trend which becomes explosive for n > 7, where fixed- 
order perturbation theory breaks down. The oscillatory behaviour is exactly 
what one would anticipate from the alternating-sign factorial growth of the 
contribution of the leading UVi renormalon, given by Eq.(26). To attempt 



15 




0.42 0.44 0.46 



Figure 2: versus as(m^). The dotted curve is the exact NNLO CORGI 
fixed-order result, and the upper sohd curve is the approximate all-orders 
CORGI resummation. 



to estimate the uncertainty in as(m^) extracted from Rr measurements we 
can use the difference between the resummed and exact NNLO fixed-order 
CORGI results to estimate the possible effects of uncalculated higher order 
terms. In Figure 2 we have plotted R^ versus as(m^). The upper solid 
curve is the all-orders CORGI result, whilst the lower dashed curve is the 
NNLO fixed-order CORGI result,. We note that the separation of the curves 
increases rapidly with increasing R^, so we are fortunate that for the ex- 
perimentally measured i?T-~0.2 the separation of the curves is reasonably 
small. Using the ALEPH data we find asim^) = 0.330lg:|]i3 from the all- 
orders CORGI result, and ^.(m^) = 0.355l[5:[Ji from NNLO fixed order 
CORGI. The corresponding results which would have been obtained by inte- 
grating up the Effective Charge (EC) beta-function for D as in Ref.|]19[| are 
a,(m2) = 0.337l|];[j}^ and ^.(m^) = 0.347l[5:o22 for the resummed and NNLO 
EC results. So, as expected, the two approaches yield similar results. If we 
evolve these ^^(m^) results through flavour thresholds up to /x = Mz using 
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the three-loop matching conditions @, 0, we find a,(M|)=0.120l[]:{j^2 
from the resummed CORGI result, and a,(M|) = 0.123io:oo2 from the 
NNLO CORGI result. Thus, we conservatively estimate an uncertainty 
Sas{Mz)^0.003. A direct plot of the resummed and NNLO results for 
versus as(M|) is given in Figure 3. The invariant mass distribution of the 
produced hadrons in r decay is well-measured experimentally ||2^, |23[| . 
We define the quantity Rt{so) as 

r (r^z/^ + hadrons; s/,ad > So) p rf_R^(s) , , 

Rr{so)= = / ds— , 37) 

1 [T^UT-eUe) JO as 

where denotes the measured inclusive hadronic spectrum. 



Rr{so) = N{\Vudf)SEw[{2x - 2a;3 + x") + ^cpRriso) + 5pc] , (38) 

with x=So/ml. The perturbative part Rr{so) can be computed from Eq.(lO) 
with the choice of weight function 

W{e) = 2x{l + e'^) - 2x\l + e^^^) + x\l - e^'^) . (39) 
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It is then straightforward to obtain contour-improved fixed-order and re- 
summed CORGI results for R{so). In Figure 4 we show the fit of the all- 
orders leading-6 CORGI resummation (solid line) to the ALEPH data for 
Rr{s) (open circles) ||2^. The resummation is fitted to the data at s = m^, 
where Rrijnl) = Rr. The CORGI coupling has a Landau pole at = Ad, 
as is apparent from Eq.(17). Fitting to the experimental value of Rr deter- 
mines Ad = 0.725 GeV, and so the resummed prediction is only defined for 
s > 0.525 GeV^. There is excellent agreement with the data. On this scale 
the fixed-order NNLO CORGI result would not be distinguishable from the 
all-orders result, and so we have not included it on the plot. 



5 Estimating the uncertainty in hadronic cor- 
rections to a(M|) 
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Figure 5: Fixed order results (starred points) for R versus different orders of 
perturbation theory at ^/s = Mt-= 1.777 GeV. The sohd hne shows R for the 
all-orders contour-improved resummation. 



In this section we wish to make use of the difference between the NNLO 
fixed-order and resummed CORGI results for R{s) in e~^e~ annihilation to 
estimate the uncertainty in a(M|), the QED coupling at the Z pole, which 
plays a crucial role in constraining the Standard Model Higgs mass from 



precision electroweak fits to radiative corrections |35|. We begin, however 



by plotting some figures, analogous to Figure 1, to indicate the performance 
of fixed-order perturbation theory versus the resummed results at various 
energies. In Figure 5 we show the all-orders CORGI leading-6 resummation 
(solid line) and fixed order results (starred points) for R{s) at ^/s = 1.777 
GeV, corresponding to rrir, so the performance can be directly compared 
to Figure 1. The only difference in the two calculations is the choice of 
weight function, W{6) in Eq.(lO). The oscillatory trend due to the leading 
ultraviolet renormalon is again evident, with wild oscillations setting in at 
n > 9 where fixed-order perturbation theory breaks down. 
In Figure 6 we present a corresponding plot at LEPl energy ^/s = Mz- 
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Figure 6: As Fig 5, but at y/s = Mz 



Clearly at the higher energy the agreement is much improved. With the fixed- 
order results exactly tracking the all-orders result for n > 4. Wild oscillations 
only set in for n > 30 at this higher energy. 

Finally, in Figure 7 we show a plot of R{s) versus ln(y^/GeV), in the range 
1 < ^/s < 91 GeV. The solid line corresponds to the all-orders resummed 
result and the dashed line to the NNLO fixed-order CORGI result. We 
assume q;s(M|) = 0.119, and evolve through fiavour thresholds using the 
three-loop matching condition [^, |3^ . 

The QED fine structure constant is extremely well-measured with 



a 



-^=a{Oy^ = 137.03599976(50) 



If one wishes to evolve a away from s = to obtain a{s) 
needs to know leptonic and hadronic corrections, 

,-1 



(40) 

however, then one 



a{s) = (1 - Aaicp(s) - Aahad(s) - Aatop{s))a ^ 



(41) 



Whilst the leptonic corrections are known at three loops and are well-determ- 



ined [36], the hadronic corrections for the contribution of the five lightest 
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Figure 7: Fixed order (dashed line) and all-orders renormalon resummations 
(solid line) for R{s) versus ln{^/s/GeV), over the range 1 < ^/s < 91 GeV. 



flavours, which we have denoted AQ;had('S), is rather poorly determined and 
has to be reconstructed from the s-dependence of R{s) using a dispersion 
relation. The contribution of the heaviest flavour A«top(s) is rather well- 
determined and can be included separately. The value of a(M|) is of par- 
ticular relevance since it limits the precision with which the unknown Higgs 
mass Mh of the Standard Model can be predicted from precision electroweak 

~^ and 



corrections Taking s 
Aatop(M|) = -0.7 
dispersion relation, 



M| we have H Aaiep(M|) = 314.98x10 



AQ;top(M|) = —0.76x10 ^. For the hadronic contribution we can use the 



A«had(M: 



J 4ml sis 



M|) 



(42) 



In Ref.[24| new exclusive data from BES-11 |^ and Novosibirsk [38| have 



been used to extract R{s) in the low energy region, with NNLO fixed or- 
der perturbative QCD used to evaluate it in the ranges 2.8 < -y/i < 3.74 
and 5 < ^/s < oo. We plan to approximate R{s) in these latter ranges 
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by the all-orders and NNLO fixed-order results for R{s), as plotted in Fig- 
ure 7. We shall use the exclusive data results as in Ref.p4|, in the re- 
maining energy ranges. Taking as{Mf) = 0.119 we shall then determine 
a{M'^) from the fixed-order CORGI results, and the all-orders leading-6 
resummed results. Since these results are contour-improved they include 
a resummation of analytical continuation terms not included in the fixed- 
order perturbative results used in 



24|. 



We are interested in establishing 
if these terms and the uncalculated higher-order corrections, as estimated 
by the leading-6 approximation, cause a significant shift in a (M|), and 
whether this has any ramifications for the constraints on Mh- In the region 
2.8 < < 3.74 we obtain Aahad(M|) = (9.5424xl0-^ 9.7075x10"^) for 
the (fixed-order, all-orders) CORGI results, and in the region 5 < y/s < oo 
we find Aahad(M|) = (170.788xl0-^ 170.635x10-^). We find correspond- 



^128.967,128.971), to be compared with 



ingly using Eq.(41), a(M|)~ 

a(M|)~^ = 128.978±0.027 quoted in Ref 0. We conclude that the analyti- 
cal continuation terms and uncalculated higher order perturbative corrections 
do not cause a significant change in a(M|), and their inclusion does nothing 



modify the conclusions of Ref [24 



6 All-orders CORGI resummations for the scalar 
correlator 

The Higgs decay width to a quark anti-quark pair will be of fundamental 
phenomenological importance. In practice the decay to a bb will be the 
dominant contribution 

m^bb) = ^MnmliMDniMl) . (43) 

Here Mh is the Higgs mass and mb{Mfj) is the running 6-quark mass. TZ{M^) 
is a coefficient function with a perturbative expansion 

7^(M^) = 1 + J2rna- , (44) 

71>0 



where the coefficients ri,r2,r3 have been exactly computed TZ can be 
straightforwardly related to the scalar correlator H5(s). One can define an 
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analogue of the vector Adler D-function so that (cf. Eq.(8)) 



Dis] 



ds 



njs) 



(45) 



This may be written in terms of the coefficient function where 

3 



D s) 



and V has the perturbative expansion 

Vis] 



1 + J2dna' 

n>0 



(46) 



(47) 



ml{s)TZ{s) can be related to m^(— s)D(— s) by analytical continuation from 
Euclidean to Minkowski, and one can write a representation of the same form 
as Eq.(lO) with 



m2(M|) 7^(M|) 



1 



de ml{e'' Ml)V{e'' Ml) 



(48) 



To proceed further we can express the running mass in terms of an RG- 
invariant mass and the mass anomalous dimension 7m (o), defined by 



(iln(m(s)) 
dhiis) 



i>0 



A+1 



(49) 



We can then write ^1 



2/ \ - 2r 

mb [s) = m^b b 



1 + ca 



l2Q 
6 



exp 



°<jx ^'^ ^^^^ + 72 - 10C2)X + ... 
+ CX) (1 + CX + C2X^ + . . .) 



470 



(50) 

The is the standard normalization of the definition of the RG-invariant 
mass rhb. One can then write a CORGI series for ml{s) ^{3) exactly equiv- 
alent to that for the moments of structure functions in Ref. ||21||, 



771^(5)1^(5) = 777^6 



470 



470 
b 



1 + cao{s) 



'1+X2al{s)+Xsal{s)+. . .+X^a^{s)+. 



(51) 
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where ao(s) denotes the CORGI couphng which is again defined in terms of 
the Lambert VT-function as in Eq.(17), and with the anomalous dimension 
present one now has = exp[(;^) + (^)](y) ''Ajjg , with d the coefficient 
di in the MS factorization and renormahzation scheme with = = s (M 
denoting the factorization scale) [^]. The exactly known CORGI invariants 
X2 and X3 follow from Eqs.(18) of Ref.pll , and allowing for the different 
definition of the anomalous dimension one needs to replace di in Ref.[pl|1 by 
. Lumping various inessential prefactors together we can define 

T{H^bb) = f^M^m^647o/ftr , (52) 
where F has the contour-improved CORGI representation, 

r = :^/ de(—^) (1 + X^a^ + Xsa^ + ^Xf )<) , (53) 
with ao = ao{e'^Mjj). 

In the scalar case one will have coefficients with the structure 

dn = 4"-iix;-i + + . . . + 4°' , (54) 

and after replacing Nf = ~ '^^)' before, one arrives at a leading-6 
term with the structure d^={—3)^~^ d]^~^^b"'~^ , with one less power of b. The 
anomalous dimension coefficients will have the structure 7^^^ = 7^"^&", but 
since the anomalous dimension 7rrt(a) does not contain renormalons there 
is no motivation for making this approximation, and it is poor in practice, 
as noted in Ref . ||25|| . Whilst an all-orders result for 7^^) does exist ^5 
shall follow Ref. ||25|| and set 7^^^ = for n > 0, retaining only 70. The all 



we 



orders result for X^^^ follow straightforwardly from d!f^\ From the large-X/ 
results of Ref.[^] for the scalar correlator one can obtain an explicit all- 
orders expression for dl^\V) (in the V-scheme) analogous to Eq.(21) in the 
vector case. For n even one has, 

4''(V0 ^ -f ^ (1 - i) C(« + DnUf-^ + (M) (i)"".*'-' 
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whilst for odd n one has, 

As in the vector case one can define a leading-6 resummation 

oo 

P(^) = 1 + E4^^(V^K, (57) 

k=l 

analogous to Eq.(22), which may be defined as a regulated Borel sum 

roc 

V^^\l/a) = 1+PV dz e-"/"[G'_(2) + G+{z)] . (58) 

JO 

Here G-{z) and Gj^{z) are the contributions to the Borel transform from UV 
and IR renormalons, respectively. One has (in the V-scheme) |2^ 



G-(-) = 



oo /' 1 \fc 



3^iA;2(l + 6z/2A;)^ 



4 1 4 °° (-1)'' 

^^^"^^ = (1 - hz/2) " (1 - hz/A) + 3£f(1-6z/2A;)^ ^^^^ 

^From these expressions one can read off the residues Aq, Bq (cf. Eq.(24)), 
and one can then calculate V^^\F)\^y and 1)^(^)1^^ using Eqs.(27),(28). 
Finally 

V^'\F)=V^^\F)\^y + V^^\F)\,^, (60) 

with F=l/ay. To perform the leading-6 CORGI resummation in Eq.(53) we 
simply need to relate ay^ and a^^^ as we did in Section 3. In the presence 
of the anomalous dimension the RS-invariant combination po iii Eq.(32) is 
replaced by the factorization scheme and RS (FRS) invariant combination 
Xi introduced in Ref. [ 



Xi = 47olnf^')-^-rfi(M), (61) 



V A 

where M is the factorization scale. Recalling that we have decided to set 
7i^^ = for z > in our leading-6 resummations , we can use Eq.(61) to 
relate Oq^^ and ay\ 

^ ^ —d['\v). (62) 



il^^ ai") 470 
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Figure 8: Fixed-order CORGI results for F {Mh = 115 GeV) in N'^LO 
perturbation theory (starred points), compared to the all-orders resummation 
(solid line) 



It then follows that the all-orders formal resummation in Eq.(53) is given 

+V^L) fl_ ^ ^d[')^v)] + (X, - ))a^ + (X3 - Xf ))a^] . (63) 
\ao 470 J 

The logarithm term arises because of the fractional power a^'^°^^. Relating 
the V^-scheme and CORGI couplings at the leading-6 level one has, 

On expanding using the binomial theorem only the terms linear in 70 are 
leading in b, the remainder should be discarded. Writing the binomial ex- 
pansion as exp[(— 470/6) In^] = 1 — (47o/6)ln5' -|- 0(70^), the result follows. 
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The same subtlety enters in deriving an analogue of Eq.(36) to generate the 
X^^) in terms of SJ^\V) explicitly given by Eqs.(55),(56). One finds 



x(^) = C 



n 



iti '\l + {A^o/b)Si'\v)aJ b \ b 

(65) 

These results can then be used to calculate all-orders and fixed-order CORGI 
predictions. We give in Figure 8 the analogue of Figs. 1,5,6 , for the Higgs 
decay width (with pre-factor set to unity) F, we set Mh = 115 GeV and 
as(M|) = 0.119. The starred points show the fixed-order CORGI results, 
and the solid line the all-orders resummation. As before the agreement of the 
highest exactly calculated n = 3 fixed-order result with the all-orders result 
is good. The fixed-order results track the resummed result up to n = 12, 
beyond which an oscillatory trend is noticeable. From Eq.(59) one can see 
that UVi and IRi renormalon singularities are present, and so the leading 
asymptotics are not dominated by UVi as in the vector case. The process 
of analytical continuation , however, serves to remove IRi 0, |] and so the 
leading asymptotics of F is expected to be dominated by the leading UVi 
renormalon, with resulting alternating-sign factorial behaviour. As discussed 



in Ref . the presence of the leading IRi renormalon in T) suggests that the 
obvious generalization of the Adler function in Eq.(45) may not be optimal, 
and an alternative is suggested. For our purposes here we are simply using 
D as a tool to compute the physical quantity F, and so this is not a problem. 



7 Analytic expressions for the CORGI contour 
improvement 

In this section we wish to point out that we can obtain explicit analytic 
expressions for the CORGI fixed-order contour improved results for R{s) 
in terms of the Lambert M^-function, eliminating the need for numerical 
Simpson's rule evaluation. From Eq.(17) we see that we can write the CORGI 
coupling clq = ao(e*^s) in terms of the Lambert ly-function as, 

^° = c[l + W{Ae^^^)] ' (^^^ 
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where 

Thus after the contour integration the Xn-i coefficient multiphes 

+ 2^/0 '^^^^^ + W^-i(A(s)e^^^)]" , (68) 

where the appropriate branches of the VT-function are to be used in the two 
regions of integration. By making the change of variable w — W{A{s)e^^^) 
we can then obtain the above integrals in the form, 

(69) 



w 



The ty-integral is elementary, and including the limits of integration, and 
noting that Wi{A{s)e~'-^'^) = [W_i{A{s)e''^'^)]* , we obtain the explicit result. 



(-IT 
AJs) = ^—^Im 



1 + W_i{A{s)eik^ ) ^^k(l + W_iA(s)e'^^)'' 

(70) 

for n > 2. For n=l we have Ai{s) = {-1/ {7iKc))Im[lnW_i{A{s)e'^^)]. We 
finally obtain the CORGI contour-improved fixed-order results in the form, 

R{s) = Ai{s) + Y,XkAk+iis) . (71) 

k=2 



In the scalar correlator case analytic results can also be obtained. The 
Xn coefficient in the CORGI series for T in Eq.(53) will multiply 



1 /-O^^ (_!)" (.-n-4^,„/h 

' ^J-^ [-Wi{Ae^^)f^°'^ [1 + Wi{Ae^^)f 

1 /-TT ( — „-n-47o/6 

271 Jo [-iy_^(/le^-^^)]^^°/Ml + W"-i(Ae^-^^)]" ■ ^ ^ 
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Here A = A{Mfj). Making the change of variable u = —W{Ae^^^) one can 
then obtain the above integrals in the form, 



'-ir 



duo 



^(-470/fc)-l 



n-1 • 



(73) 



These integrals may be evaluated in terms of the Hypergeometric function 
F{a,b; c; z) Inserting the limits of integration we obtain the explicit 

result for > 0, 



An{M] 



F n — 1, 



Im 
1 



l2a 

b 



b b 

We can finally write the CORGI contour-improved result in the form 



(74) 



r = Ao(Mi) + Y^XkAkiMfj) , 

k=2 



(75) 



with 



Ao(M 



Im 



-^{W.i{Ae'''n~^'"' + 
470 



(76) 



8 Conclusions 

In this paper we have focussed on obtaining exact fixed-order, and leading-6 
estimated all-orders results for various inclusive QCD Minkowski observables 
, related to the vector correlator. These could be expressed as a contour inte- 
gral of the suitably weighted Euclidean Adler D(—s) -function in the complex 
energy squared plane. -D(s) is truncated at some fixed-order and the inte- 
gral performed numerically as described in Section 2. In this way contour- 
improved predictions are obtained, in which an infinite subset of known and 
potentially large analytical continuation terms are resummed to all-orders. 
By employing the CORGI approach, as discussed in Section 3, we could 
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further resum to all-orders the complete set of ultraviolet logarithms involv- 
ing s, which build the s-dependence of -D(s), avoiding any dependence on 
an arbitrary renormalization scale /i. The remaining approximation is the 
missing higher-order CORGI invariants Xi {i > 2), which remain unknown 
since the perturbative coefficients have only been calculated to NNLO so 
far. We used the so-called leading-6 approximation to estimate these. Us- 
ing the exact large- A^/ all-orders results for D{s) |0, H we were able to sum 
the CORGI series to all-orders in terms of a sum of exponential integral 
functions, corresponding to the contributions of the ultraviolet and infrared 
renormalons in the Borel plane. This was technically far more straightforward 
then previous analogous resummations of the Effective Charge beta-function 
for D{s) |TB|, |T3. By comparing the NNLO fixed-order CORGI results to 



the all-orders resummations we estimated the uncertainty in as{M^), ex- 
tracted from experimental measurements of Rr, to be 5a;s(M|)^0.003. We 
also showed that using all-orders and fixed-order contour-improved CORGI 
results for R{s) gave results for the hadronic corrections to the QED coupling 
a(M|) which did not differ significantly from those obtained using conven- 



tional fixed-order perturbation theory in Ref.[^. Finally we showed how 
recent all-orders large-iV/ results for the scalar correlator could be used to 
perform analogous resummations for the Higgs decay width to a heavy quark 
pair. We finally noted that the CORGI contour-improvement for the i?-ratio 
can be written in analytic form in terms of the Lambert VT-function, and 
for the Higgs width in terms of hypergeometric functions and the Lambert 
VT-function, thus avoiding the need to use a numerical Simpson's Rule eval- 
uation. 



There are many possible investigations to be pursued using these meth- 
ods. In particular it would be interesting to investigate hadronic corrections 
to a(M|) using the resummed results for R{s) in lower energy ranges where 
conventionally inclusive or exclusive data has been used. The scalar correla- 
tor results could also be used to investigate more carefully the uncertainties 
in estimates of the bottom quark mass using Sum Rule techniques. 
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